Spinor fields in Causal Set Theory 
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Abstract 



The goal of this paper is to define fermionic fields on causal set. This is done by 
the use of holonomies to define vierbines, and then defining spinor fields by taking 
advantage of the leftover degrees of freedom of holonomies plus additional scalar fields. 
Grassmann nature is being enforced by allowing measure to take both positive and 
$h | negative values, and also by introducing a vector space to have both commutting dot 

product and anticommutting wedge product. 



>■ ' 1. Introduction 

This paper is the fourth paper in a series of papers whose purpose is to define quantum 
. field theory in a causal set. A causal set (for recent review see Ref [5] ) is a way of discritizing 
spacetime while retaining strict relativistic invariance. This is done by removing coordinate 
system, and viewing spacetime as a partially ordered set, where partial ordering is physically 
viewed as light cone causal relation between different elements of the set, that are viewed as 
events. Discritization is enforced by the postulate that any interval, a(p, q) = {r\p -< r -< q} 
contains only finitely many points. One of the challenges of the theory is the fact that fields, 
as they are normally denned, have Lorentzian distances, and so do the derivatives. Hence, 
in order to be able to do physics in causal set framework, one has to rewrite physics in a way 
that makes no reference to coordinates. Scalar fields remain to be scalar valued function, 
vector fields are redefined in terms of holonomies, while gravitational field is redefined in 
terms of causal relations themselves. In papers [3], [I] and [2] a proposed way of estimating 
Lagrangian densities of these fields were given. The goal of this paper is to do the same for 
fermionic field. 

This paper is based on "translating" the model of Ref pQ from manifold to causal set 
context. The goal of Ref pQ was to re-define fermions in a way that adresses the following 
two issues: 

1) Grassmann nature of the field 2) The presence of vierbines and whether or not they 
should be viewed as fields 
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According to this approach, we substitute vierbines with a set of four vector fields that 
are no longer assumed to be orthogonal to each other. Extra degrees of freedom that are 
gained by that are re-interpretted as fermionic field. This implies that, indeed, vierbines 
are fields since both vierbines as well as spinor fields are different degrees of freedom of the 
same thing, so if the latter is a field, the former is a field as well. Furthermore, Grassmann 
nature of fermionic field is enforced by equipping that space with both commutting dot and 
anticommuting wedge product, together with cleverly defined measure that has both positive 
and negative values. 

The goal of this paper is to rewrite the results of [1] while replacing the four vector 
fields with four holonomies. The challenge is to be able to express the contractions of vector 
fields itself as well as their derivatives in a coordinate independent way, based on holonomies 
alone. 

2. Quick review of [lj 

Since this paper is based on results of p] it is best to summarize these results before 
proceeding with this paper. 

In Ref [T] fermions were described in terms of the four orthonormal vector fields, for the 
case of standard, coordinate-based manifold. It was done based on the following three steps: 



PART 1: DEFINE GRASSMANN NUMBERS AS LITERAL MATHEMATICAL OB- 
JECTS 



We introduced a space equipped with commutting • product and anticommutting A 
product, together with measure £ that has both positive and negative values, in such a way 
that J (d6i A ... A d9 n ) ■ (6 ai A ... A 6 am ) obeys expected properties of Grassmann integral. This 
allows to define Grassmann numbers in a literal sense. 

This can be done by enforcing the following two restrictions on measure £ : 




(1) 




(2) 



For the purposes of the other two steps, it is convenient to define 



dS(x) 



(3) 



dx 



which can be altered, for the purposes of rigour, to 



a 



1.5 



ax) 



.0.5 



e 



—ax 



(4) 



7T' 



Expected properties of dot and wedge product are the following: 
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X • X — 1 

y ■ ( x A y) = -y ■ (y A x) = x 



(5) 
(6) 



(y A z) ■ (x A y A 5) = x (7) 
z-(iAj/Az) = ^Az. (8) 

This, of course, is generalized to the products of larger number of multiples. 

PART 2: COME UP WITH MORE GEOMETRIC DEFINITION OF FERMIONS IN 
A TOY MODEL WHERE THERE ARE NO GRASSMANN NUMBERS. 

Introduce fields A^, B^, C M and D^. Define ftf through fg to be unit vectors in respective 
directions of N 1 through D^: 



A» £?m 

= p 1 = f» = f» = (a) 

h v^a: 1 2 3 w 

Use Gramm Schmidt process to define a set of orthonormal vectors at each point, e^, 
el, e%, and eg that are interpretted as vierbeins: 



^ - ^mt (10) 



4(A, B) = B e °^ e ° (11) 

^{BP-e%B a e p Q ){B p -e P B p e 0p ) 



S(A BC)= Z = 6 °^ e ° = ^ (12) 

V(<> - e?C 7 eg - e\C s e\){C p - e» C 5 e 0p - efC e e lp ) 



— p a D — Dvp^ — p 1 D p^ 
e£(A,B,C,D) = "° 1/3 1 2 7 2 (13) 

yj (Dp - e s D s e p - e\D e e{ - 4d^){Dp - e*D x e p - e\D n e{ - e%D^ 2 ) 

Introduce two scalar fields \p an d Xa, and then define fermionic field to be a rotation of 
Xp \ u i) + Xa \ v i) from the /^-based to e^-based reference frame: 



°' \a\* \b\' \cy \D\' v 1 

= (ex V {-\(\n(e- l (A,B 1 C 1 D)f(A 1 B,C,D)))^}) ij ( Xp 5{+XaSi) . 

In other words, while are viewed as vierbines, extra degrees of freedom given by f£ 
allow us to define fermionic field. At the same time, remembering that we have obtained 
from fg by Gramm Schmidt process, the fermionic field is viewed as function of f£ alone: 

- /rw, f%(x), fg(x)) 

PART 3: COMBINE PARTS 1 AND 2 TO OBTAIN LITERAL DEFINITION OF 
FERMIONIC FIELDS. 

Introduce unit vectors f\ through r 8 remembering that we would like r 1; f 3 , f 5 and r 7 
to stand for real parts of tpi through ip^ while f-i, r 6 and f$ to stand for imaginary parts 
of ^1 through ^4 Based on this, define the field 



tfiXp, Xa, A, B, C, D) = nRe^ixp, Xa, A, B, C,D)) + if 2 Im(ij) 1 (Xp, Xa, A B , C, D)) 
+f 3 Re(i/j 2 (Xpi Xa, A, B, C, D)) + ifJm(ip 2 (x P , Xa, A, B, C, D)) 
+f 5 Re(i/j 3 (x P , Xa, A, B, C, D)) + ir G Im(ip 3 (x P , Xa, A, B, C, D)) 

+f 7 Re(Mx P , Xa, A, B, C, D)) + if 8 Im(Mx P , Xa, A, B, C, D)) (15) 

Define a measure A(/o t , /f , f£ , f 3 ) in such a way that its "derivative" with respect to 
transformation (/^(x), fi(x), /^(x), f 3 (x)) — > ^ would give us an expected measure £(^1) 
with both positive and negative values as needed for the definition of Grassmann integral 
(see part 1) : 



/ A B C D 

KXp, Xa, A, B, C, D) = £ $^ XptXa — — — — 1 x 

r / A 1 B' C D' 

x lim^oefi- 1 jxp, Xa, A ', B' , C , D' \ip(x' P , X a , rjry 7^77 , t^jt , 7^7; 

, / A 5 C D \ , 

^l^' Xa 'R'p'M'MJl <eAVa 



A' 5' C" D' 



A B' C D' 



^afj, \ 



I AT 1ST ICT \D' 



-a/i 



a V|A'|' |C"|' \D'\ 
A B' C D' 



- e: 



A B C D 

\Ay\Wy]cy\D\ 



I AT \B'V ICT ID'I 



< e 



Here \i is a usual Eucledian measure on M 18 

We then use this information to rewrite f[di/)](...) in terms of f[dfo][dfi][df2][df £](... ) 

The bottom line is that in Ref [1] fermionic field theory was re-done in terms of four vector 
fields A 1 , B^ 1 , and D^. The goal of this paper is to "translate" this theory into causal set 
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context by replacing these four fields with holonomies, thus removing the Lorentzian index. 
As with other papers (see [3J, [3] and [2]) the ultimate challenge is to express the Lagrangian, 
which involves contractions of derivatives, in coordinate-independent way. 

2. Computting key terms of the Lagrangian 

2.1 General techniques 

Since in our philosophy there are no such things as coordinates, we will interpret vierbines 
eo, ex, e<i and as four separate fields. Our Lagrangian is 

£ = ?W + TV* ^ («(«V fc „ - cU M ) + e£ el d a e™) . (18) 

All terms of the above Lagrangian can be presented in the following three forms: A^B^ , 
A^dptfi, and A' 1 B u d u C ll . Thus, our goal is to present each of them as a function of holonomies 
corresponding to each of these fields: 



(A»B,)(p)ttf(a,b;p) (19) 
A^<t>&g{a,(j)- V ) (20) 

(A»B»d u C,)(p)^h(a,b,c;p) (21) 

We will use the same technique we have used in Ref [2],[3],[3] for bosons. In these 
papers, a region called Alexandroff set, a(p,q) = {r\p -< r -< q}, was used in finding the 
approximation to Lagrangian density. The basic steps were as follows: 

a) For any arbitrary Alexandroff set a(p,q), define £(F,a(p,q)) to be a proposed ap- 
proximation of Lagrangian density on the interior of that Alexandroff set. 

b) Establish a criteria of "reliability" of each of the Alexandroff sets by finding a way 
to measure "fluctuations" of the characteristic behavior of a given field in their interio, 
Fluc(F, a(j>, q)), and the bigger the fluctuations are the less reliable Alexandroff set is. 

c) For any point p define a point q(F,p) in such a way that it minimizes Fluc(F, a(p, q)) 
under the constraint that a(p, q) has more than certain fixed number of points (the latter is 
needed in order for statistics to be reliable) 

While the above basic steps are the same, there are two very different directions of 
implimenting them. 

l)Integration approach. This apporach was used in [3j and [3]. According to this ap- 
proach, part a is carried out in such a way that C(F,a(p,q)) will give exact value of La- 
grangian density, as long as spacetime is flat and all fields are linear. While this in itself will 
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be obvious advantage as compared to the next approach, the disadvantage is that, as seen 
in [3] and [I], in order to enforce such a criteria, we need to define Lagrangian to be a linear 
combination of two or more terms with non-integer coefficients that are very different from 
any numbers otherwise recognized. This adjustment is needed in order to make sure that, in 
locally flat coordinates, we have a = — b in C = oTqq + b(Tu + T22 + T33) This immediately 
raises a question as to how come we have to enforce relativity "by hand" if causal set, being 
based on -< alone, should be inherently relativistic. 

2)Maximization/Minimization approach. This is an approach that was used in [2] Ac- 
cording to this approach, in part a we don't bother making sure that the Lagrangian density 
we propose actually matches the one expected. Instead, we just take one relevent term and 
"forget" about any others. What we will find is that after having done part c, the other terms 
that we dropped will no longer be an issue because we have selected q in such a way that 
Alexandroff set a(p, q) will pick out a frame in which the omitted terms are 0. For example, 
in case of scalar field, we can focus on its time derivative; in part c we will select a frame in 
which t axis points in the direction of gradient of scalar field, thus time derivative will be 
all we need for the Lagrangian density. The advantage of this approach is that we no longer 
have to adjust coefficients since we only have one term to deal with. The disadvantage is 
that this makes Lagrangian density dependent on the maximization/minimization procedure 
as opposed to being universally defined for arbitrary Alexandroff set. 

We will now use both approaches to compute /, g and h. 
2.2 Maximization/minimization approach for fermions 



I will introduce holonomies in a way similar to paper |4j. Namely, I will associate the 
above vector fields with scalar fields of pairs of points that are defined in the following way: 

a(r, s) = (a" - r") (^(0) + ± (r" + s v ) d u A^) + 0(r 3 ) 
b(r, s) = K - r") (£„(0) + \ if + s") d v B^) + 0(r 3 ) 
c(r, s) = {s» - r») (C,(0) + \ (r" + s v ) 3„C M ) + 0(r 3 ) 
d(r, s) = (*" - r") (£>„(0) + £(/" + *") d u D^ + 0(r 3 ) . (22) 

The generic things I would like to compute for causal sets is a term of the form A tl B u d l ,C^ 
as well as B^. 

Let us first compute gauge latter terms. If we fix point p and think of a(p, r) and b(p, r) 
as scalar functions of r alone, then we can use the result of the scalar field section and say 



\A\ 



\B\ 



da(p, r) 




Q r ii 


= ( 


db(p, r) 









k A 1/d ■ 1 ( ^ 

— I mm max \a(r,s) 

Vq J V(a(p,q))>Vo r,s£a(p,q) 
l/d 

min max \b(r,s) 



kd 
Vo 



\A + B\ 



d(a(p, r) + b(p, r)) 



Q r fj, 



V(a(p,q))>Vo r,s£a(p,q) 
kA 1/d 

Vo 



mm max 

V(a(p,q))>V r,s£a(p,q) 



(23) 
(24) 

\a(r,s) +b(r,s)\ . (25) 
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This tells us that 



A»B li = \{\A + B\ 2 - \A\ 2 - \B\ 2 (26) 
(k \ X ^ d 

— h\ — ) f min max \a(r, s) + b(r, s)\ 

2 \V ) \V(a(p,q))>Vor,sea(p,q) 

min max |a(r, s)\ — min max \b(r, s)\) . (27) 

V(a(p,q))>Vo r,s£a(p,q) V(a(p,q))>Vo r,s£a(p,q) J 

Due to the fact that we seen in scalar section that the same equation applies both to spacelike 
and timelike gradients of 4>, it is clear that the above question would be the same if we were 
replace A 1 with one of the spacelike vectors. Thus, 

5" C„ = i( — ( min max \b(r, s) + c(r, s)\ (28) 

M 2 V V J V V(a(p,q))>V r,sea(p,q) 

— min max \b(r, s)\ — min max |c(r, s)\) . 

V(a(p,g))>Vb r,sea(p,q) V(a(p,q))>V r,s£a(p,q) J 

Now, this can also give us the value of A 1 <9 M Namely, I should just replace in my mind 
b(r, s) with <p(s) — 0(r), which means replacing with d^<p. Thus, we have 

A " d ^ = \\4) ( w mi " max |a(r, S )+0( S )-0(r)| (29) 

V v J v V(a(p,q))>Vo r,s£a(p,q) 

— min max |a(r, s)\ — min max |0(s) — <f>(r)\ ) . 

V(a(p,g))>V r,s€a(p,q) V(a(p,q))>V r,s£a(p,q) / 

Again, the same equation applies if we replace A* 1 with something spacelike: 

^^ = Mt7 Lr* max |6(r, S ) + 0(,)-0(r)| (30) 

— min max \b(r, s)\ — min max |0(s) — <f>(r)\ ) . 

y(a(p,g))>V() r,s£a(p,q) V(a(p,q))>Vo r,s€a(p,q) J 

Now let us try to get A* B v d u C^. 

Suppose we would like to maximize the fluctuation of a(r, s)b(r, s)c(r, s)d(r, s) Ideally, if 
we could minimize the fluctuation of each of individual four multiples, it would automatically 
minimize the fluctuation of the product. I claim that we can do just that! Based on the results 
of scalar part, if gradient of (p is timelike, in order for fluctuations of </> to be minimized it has 
to be parallel to the axis of Alexandroff set; on the other hand, if gradient is spacelike then 
in order for fluctuations of 4> to be minimized, gradient should lie on equator of Alexandroff 
set. This means that in order for fluctuations of a(r, s) to be minimized, A^ should point 
along the axis of Alexandroff set, while in order for fluctuations of b(r, s), c(r, s) and d(r, s) 
to be minimized, then B^, and _D M should lie on equator of Alexandroff set. Now, the 
amaising thing is that orthogonality condition tells us that, not only these four statements 
are compatible, but in fact if the condition about the timelike vector A 1 is met, it forces 
the conditions about the three spacelike vectors and to be met as well! This 

means that if we minimize the fluctuations of a product a(r, s) b(r, s) c(r, s) d(r, s) we would 
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also minimize the fluctuations of each of the four multiples individually, which means that 
we would completely specify the paramenters of my Alexandroff set. 

Now suppose we would like to select points r±, r 2 , r 3 and r 4 in such a way that minimizes 
fluctuations of a(r\,r 3 ) b(r 2 , r 4 ). Again, we can separately maximize each of these two multi- 
plets. To maximize a(n, r 3 ) we have to set r\ = p = (— \ r, 0, 0, 0) and r 3 = q = (| r, 0, 0, 0). 
In order to maximize 6(r 2 ,r 4 ) we have to select r 2 = (0, |r, 0,0) and r 4 = (0, — \ t, 0, 0, 0) 
These two conditions are comparable with each other. Now suppose that we instead decided 
to maximize fluctuations of b(ri,r 3 ) c(r 2 ,r 4 ) Again, we can simultaneously maximize fluctu- 
ations of each of these fields. This time, since both of these fields have spacelike gradient, 
each one will be maximized by selecting points on the equator. In particular, to maximize 
K r i) r 3) we se t r i — (0, — \ r(p, q), 0, 0) and r 3 = (0, | r(p, q), 0, 0) and in order to maximize 
c(r 2 ,r 4 ) we set r 2 = (0, 0, \ r(p, q), 0) and r 4 = (0, 0, — \ t(p, g), 0) Again, we get the same 
square look as we had in gauge case. We notice that the square loop for the case of one 
timelike and one spacelike field is the same as the square loop for two spacelike fields, except 
that t axis was replaced by the y axis. So, we will save ourself time and say that in order 
to maximize u(ri, r 3 ) v(r 2 , r 4 ) we have to set ri M = — | r, r 3fM — \ r, V^r 2 ^ = \ r, 
V^r 4M = -\t regardless of whether the fields are spacelike or timelike. 

Now if w is a third holonomy, then based on our results from gauge part, remembering 
that all vectors are unit vectors, we have 

w(n, r 2 ) + w(r 2 , r 3 ) + w(r 3 , r 4 ) + w(r 4 , n) = r 2 (d p W a - d a W p ) . (31) 

where p and a are directions corresponding to U and V. Remembering that U and V are of 
unit length, this can be re-written as 

w(n,r 2 ) + w(r 2 , r 3 ) + w(r 3 , r 4 ) + w(u, n) = r 2 U» V u (d^W v - d u W^). (32) 

Now since = we know that 

U» V u d v W„ = -W V u d u U^ ■ (33) 
Thus the above expression becomes 

w(r u r 2 ) + w(r 2 , r 3 ) + w(r 3 , r 4 ) + w(u, n) = r 2 (p, q) (U» V u d^W u + W V v d v U„). (34) 

In order to remember what we were minimizing or maximizing, this can be written in a 
more complete form as follows (here we used V(a(p, q)) = k d r d (p, q)): 

{U» VdnWv + W V v d v U^) 

= (— ) 1/d min max (w(ri, r 2 ) + w(r 2 , r 3 ) + w(r 3 , r 4 ) + w(r 4 , n)) (35) 

Vq V(a(p,q))>Vo ri,r2,r3,r4€a(p,q) 
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We can now permute this equation to get another two equations: 



= (— ) 1/d min max (u(ri, r 2 ) + u(r 2 , r 3 ) + u(r 3 , r 4 ) + w(r 4 , n)) (36) 

Wq VT(a(p,q))>Wo ri,r2,r-3,r4Ga(p,g) 



and 



tf" 9 M K + V t/" d„wg (37) 
= (— ) 1/d min max (v(r u r 2 ) + v(r 2 , r 3 ) + v(r 3 , r 4 ) + v(r 4 , n)). 

Uo U(a(p,q))>U ri,r 2 ,r3,r 4 ea(p,g) 

By subtracting the third equation from the sum of first two equations, and then dividing 
the whole thing by 2, we get 

V v d^Wy = \ ((^) 1/d min max (w(r 1 ,r 2 ) + w(r 2 ,r 3 ) + w(r 3 ,r 4 ) + w(r 4 ,r 1 )) 

Vq V(a(p,q))>V r 1 ,r 2 ,r :i ,r 4 Ga(p,q) 

+ {—) 1/d min max (u(ri, r 2 ) + u(r 2 , r 3 ) + u(r 3 , r 4 ) + u(r 4 , n)) 

Wo W(a(p,q))>W ri,r 2 ,r 3 ,r4Ga(p,q) 

-(— ) 1/d mm max Mn, r 2 ) + v(r 2 , r 3 ) + v(r 3 , r 4 ) + v(r 4 , n))), (38) 

Uq U(a(p,q))>Uo ri,r2,r3,r4ea(p,g) 

which is our desired term. 

2.3 Integration approach for fermions 

Lets start with U^V^ terms. 

Suppose we have two fields XJ^ and given by two holonomies u(r, s) and v (r, s) We 
select an Alexandrov set given determined by poins p -< q in which the holonomies are 
assumed to be linear. As with other fields, we first pretend that we DO have coordinates, 
which we treat as uknowns, to get coordiniate-free estimate for Lagrangian density. We 
will then formally use that Lagrangian density as a definition of Lagrangian density of all 
Alexandrov sets, including non-manifoldlike. 

We will set time axis to go from p to q. 

_ u(p,g) _ y(p, q) 

Uo — —, r , vq - —, r • V>y) 

T[p,q) T[p,q) 



For simplicity, we will just indicate r(p, q) = t. 

f d d ru(p,r)v(p,r)= [ (r" r» U, V v ~ U V )d d r = r d+2 (U V (I d0 -^)+U k V k I dl ) 



(40) 
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This gives us 



U * V k = ^-i& I d d ru(p,r)v(p,r)-U V (I d0 - ^)) (41) 

*dl T J a(p,q) ^ 
If 1 k d 

d d ru(p,r)v(p,r) u(p, q) v(p, q) (7 d0 - ~r)) ■ 



hi r d + 2 J aM t 2 
Thus, we get 

U» V, = 1 u(p, q) v(p, q) - U k V k (42) 
= J- (\(IdL + ^-Ido)u(p,q)v(p,q)--j^ [ d d ru(p,r)v(p,r) 

*dl \T ^ T Ja(p,q) 

Furthermore, we can use the expression for U^V^ to get the terms involving derivatives 
of scalar fields in a Lagrangian by replacing v(r, s) with <f>{r) — 4>{s): 

U" <9 M (43) 

= (hi + y - ho) u(p, q) (<P(q) - 4>{p)) - / d d r u(p, r) (0(r) - <f>(p)) 

1 dl \i ^ i J a(p,q) 

Now let us move to the more difficult issue of computting e^e^d v e CfM terms. 
Our plan is the following: 

PART 1: By using orthonormality of e k , show that e a ^ e b d v e cll can be expressed as 
linear combination of the terms of the form Ei mn = ef" e m u {d^e nu — d v e njl ). This would 
simplify the situation tremendously since the latter somewhat resembles gauge theory which 
we already know how to do. 

PART 2: Find a way of computting E a j, c in a coordinate-free setting of causal set. Even 
though, as remarked above, the resemblence to guage theory should make it easy, there are 
still difference with gauge theory, including the fact that we have 3 holonomies rather than 
1, which makes it somewhat difficult. But orthonormality of these three holonomies will 
definitely help us pass through this. 

PART 1 

If we expand the above expression and switch the dummy indices fi and v on the second 
term we get 

E a bc &b d^C cy 6 a 6^ 0^6 CL , . (44) 

If we now apply the fact that 9 M (e a M e CM ) = d^rjab = to the second term of above equation, 
we get 

Kbc = e b u e/ d^e cv + e c " e^d^ . (45) 
By permuting the indices, we get 

E bca = tc e 6 " d^e au + e/ e/d^e bu (46) 
E cab = e a u e/ + e b u e/^e c „ . (47) 
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From these expressions it is easy to see that 

E abc + E bca - E cab = 2 e/ e b » d^ a v . (48) 
By switching a and c and dividing the expression by 2 we get 

e/ e b » d^ v c = X - (E cba + E bac - E acb ) . (49) 

Thus, the problem of computing e a v e b ^ d^e c u reduces to the problem of computing E cba for 
causal set, as desired. 

PART 2 

Let us first make a heuristic argument to make the best guess as to what kind of integral for 
us to try. 

We know from the case of electrodynamics that 

e c (p, r) + e c (r, s) + e c {s, p) = (r" - jf) {s u - f) (<9 M e c , - d u e Cfl ) . (50) 

Furthermore, if we multiply the above expression by e a (p, r) e b (p, s), we would get a pro- 
jection onto fi = a and v = b thus singling out a M b u (<9 M e CJ y — d v e cp ) The reason we can do 
this easily is because we already know the coordinates of p which leaves us with only one 
unknown coordinate. Of course, the above is just a heuristic argument so now we have to 
perform in order to see that it works. So we would like to compute 

/ d d rd d s(e c (p,r) + e c (r,s) + e c (s,p))e a (p,r)e b (p,s) . (51) 

Ja(p,q) 

Let us denote <9 M e C j, — d u e c ^ by H^. Thus, according to the current notation, 

E abc = aWH^ . (52) 

Then, after we substitute 

e c (p, r) + e c (r, s) + e c (s,p) = (r» - p") (r" - p v ) H^ c , (53) 
our expression becomes 

d d r d d s H^ c (r" - P n (s» - p v ) (r a - p)(s b - p b ) . (54) 



/ 

Ja( 



>a(p,q) 

If we now expand the parenthesis while dropping all the odd terms, we will get 

/ d d r d d s (H^ c r» s u r a s b + H Cflu r» p u r a p b + E Cflv jf s u p a s b + H^ c p» p u p a p b ) . (55) 

Ja(p,q) 

In the first term, in order for the integration over r to be non-zero, we should have fx = a. 
Alson, in order for the integration over s to be non-zero, we should have v — b. Since a and 
b are known while fx and v are dummy indices, the first term becomes H cab (r a ) 2 (s 6 ) 2 . Now, 
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by substituting p u = — \b v o and p b = — | S b o, the second term becomes ^ H^ 0c 5 b or^ r a . 
Again, in order for the integration over r to be non-zero we have to set /i = a which means 
that the second term becomes ^ H a0c 5 b o (r a ) 2 . By the similar argument, the third term will 

become ^ H obc 5 a o(r b ) 2 . Finally, p M and p v in the last term will set n = v = which would 
mean that H^ uc = H 00c = by antisymmetry. So the last term drops out. Thus, the final 
expression becomes 

/ d d r d d s (E abc {r a f {s b f + ^ 5 b E a0c (r a ) 2 + ^5 a E 0bc (s b ) 2 ) . (56) 

Ja(p,q) V 4 4 / 

Now from dimensional analysis we know the following: 

d d r = k d T d (57) 

a(p,q) 



I 

J a 



d d r (r») 2 = I d ,r d+2 , (58) 



where by cylindrical symmetry I a = I\ whenever a ^ 0. 

Substituting these we get the final answer for our integral: 



/ d d rd d s(e c {p,r) + e c (r,s) + e c {s,p))e a {p,r)e b {p,s) (59) 

Ja(p,q) 

— T 2d+A (^E abc I a h+ EaOc la $°b + ~^ E 0bc I b 5° a j . 

Now, in the situation where both a and b are non-zero (hence denoted by i and j) and 
remembering that by cylindrical symmetry (thus saying that I k = h for all k > 1 I get 

J d d rd d s (e c (p, r) + e c (r, s) + e c (s,p)) e^p, r) ej (p, s) = r 2d+4 E ijc h 2 . (60) 

Ja(p,q) 

This gives me 

Eijc = T 2 2d+ 4 / d d rd d s( y e c ( y p,r) + e c ( y r,s) + e c ( y s,p))e i ( y p,r)e j ( y p,s) . (61) 

h T Ja{p,q) 

On the other hand, if we consider a = and b > (and to stress that 6 > we will denote 
b by &) we get 

( d d r d d s (e c (p, r) + e c (r, s) + e c (s, p)) e (p, r) e k (p, s) = T 2d+A E okc h ( J + ^) . (62) 
Thus, we will get 

E ok c = ; j- fc— / d d rd d s(e c (p,r) + e c (r,s) + e c (s,p))e (p,r)e fe (p,s) . (63) 

r 2 ^ 4 /! (7 + f ) 
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Thus, the combination of results from parts 1 and 2 tells us that 

1 

2 

where 



e a v e b » d^ c = - (E cba + E bac - E acb ) , (64) 



E ijc = j 2 \ d+4 / d d rd d s(e c {p,r) + e c (r,s) + e c {s,p))ei{p,r)ej{p,s) , (65) 
when both i and j are non-zero, and 

E okc = -E k0c = — 1 j— / d d rd d s{e c {p,r) + e c (r,s) + e c (s,p))e {p,r)e k (p,s) . 

T h \h + 4 J Ja(p,q) 

(66) 

This determines e a y ef, M 9^e^. 
2.4 Definition of f, g, and h from both approaches 



Before we proceed to the next section, it is best to summarize the results of the previous 
two subsections. Since in the causal set there are no coordinates, it is not correct to say that 
we have found an approximation to the contractions of derivatives, since the latter is not 
defined. Rather, we have found functions f(u,v;p), g(u,v;p) and h(u,v;p) that happened 
to have a feature that in a special case where the causal set happens to be manifoldlike, the 
following approximations happen to hold: 



(A»B,)(p) » fi{a,b;p) » f m (a,b;p) (67) 

A^d^ w gi(a,(p;p) w g m (a,<p;p) (68) 

{A> i B v d v C li ){p) w h t (a,b,c;p) w h m (a,b,c;p) (69) 

We will now rewrite down the definitions of above /, <? and /i which would be the key 
results of the previous sections, and will be used as a part of definition of Lagrangian.The 
/, g and h that are based on maximization/minimization approach will be denoted by f m , 
g m and h m . The /, g and h that are based on integration approach will be denoted by /j, 
gi and h^. Here, p will be a point of reference, while q(u, v; p) will be a point that would 
minimize fluctuations of u and v on a(p, q) 



fk d \ 1/d / 

f m {u,v-p) = \{ — \ I min max \u(r, s) + v(r, s)\ 

\ V J \ V(a{p,q(u,v;p)))>V r,s£a(p,q(u,v;p)) 

— min max \b(r, s)\ — min max 

V(a(p,q(u,v;p)))>Vo r,s£a(p,q(u,v;p)) V(a(p,q(u,v;p)))>Vo r,s£a(p,q(u,v;p)) 
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fk d \ 1/d ( 

gm{v,<t>) = \[^) ( min max \v(r,s) + <f>(s) - <j>(r)\ (71) 

\V J V V(a(p,q))>V r,s&a(p,q) 

— min max \v(r, s)\ — min max \(f>(s) — <j>(r)\ ) . 

V(a(p,q))>Vo r,s€a(p,q) V(a(p,q))>Vo r,s£a(p,q) J 



h m (u, v,w) = \ ((— ) 1/d _ min _ max / s (tu(ri, r 2 ) + tu(r 2 , r 3 ) + w(r 3 , r 4 ) + w(r 4 , ri)) 



Vo V>(p,q))>V 



+ (— min max (w(ri, r 2 ) + u(r 2 , r 3 ) + w(r 3 , r 4 ) + w(r 4 , ri)) 

KKo W(a(p,q))>W ri,r 2 ,r3,r 4 ea(p,q) 

-(—) 1/d min max (v (n, r 2 ) + u(r 2 , r 3 ) + v(r 3 , r 4 ) + v(r 4 , n))), (72) 

Cq ^(a(Pi9))>^0 ri,r2,r 3 ,r 4 £a(p,q) 



1 / 1 , , fc, 



/,(>. :■:/•; = — ( (7 d i + ^ - I dv } ti{l).(jl n. r.j))) r(/;. </(//. r: //)) 

d d rn(p,r)f(p,r) ) (73) 



1 



r d+2 



a(p,g(«,u;p)) 



(v, 4>\p) = (J^ (hi + -£- ho) u(p, q(<f>; p)) (0(g(0; p)) - 0(p)) 

/ d d ru(p,r)(0(r)-0(p))) (74) 

hi(a, b, c; p) = ^ (£(c, 6, a; p) + a, c; p) - E(a, c, b; p)) , (75) 



where 



E(b, c, d- p) = \ d+A J d d r d d s (d(p, r) + d(r, s) + d(s, p)) b(p, r) c(p, s) , (76) 

1 l r J a(p,q(b,c,d;p)) 

when f(b, b) < 0, /(c, c) < and f(d, d) < (which is a causal set way of saying that all 
three holonomies correspond to spacelike vectors) 



and 



E(a, b, c) = -E(b, a, c) = ) j— [ d d r d d s (c(p, r)+c(r, s)+c(s, p)) a(p, r) 6(p, s) 

T 2d + 4 /j. (7 + -f) Ja(p,q) 

(77) 
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if /(a, a) > 0, while f(b, b) < and /(c, c) < (which is a causal set way of saying that 
a corresponds to a timelike vector while b and c correspond to spacelike vectors). 



3. Causal set verion of Ref [TJ 

In the previous section, we have found a way to write A^B^, A^d^cj) and A^ B h 'd v C M as 
functions of holonomies corresponding to these vector fields, 



(A^ip) « /i(o,6;p) « / ro (a,6;p) (78) 
A^d^m gfafap) » g m {a,(f);p) (79) 

{A^B^dvC^ip) « Mo,6,c;p) « h m (a,b,c;p) (80) 

where indeces z and m correspond to maximization/minimization approach while index 
i corresponds to integration approach. From now on we will drop these indeces and just use 
/, g and /i. Reader can insert either of these indeces depending on their preference of which 
approach to use. 

We will now use /, g and h to write down a complete fermionic Lagrangian. 

In order to make notation more intuitive, we will define dot products of holonomies 
(which is not to be confused with a dot product used on Grassmann space) as follows: 

(a-b)(p) = f(a,b;p)) (81) 
((a-b)c)(p,q) = (a-b)(p)c(p,q) (82) 
Thus, in a special case of a manifold, 



and 



((a-b)c)(p,q) = {A»B,){p)C u {q» - p») = c(p, q) (84) 

Now, in the same way as A^ through were just described as holonomies, we will also 
introduce a modified Gramm Schmidt process that would define orthogonal basis in terms 
of holonomies as well: 
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(eoM))(r, fl )= ^) (85) 
V( a ' a )( r ) 

/ / ny \ b(r,s)-eo(r,s)(eo-b))(r) 

[eAa, b))(r, s) = — r = (86) 

>/((6-(eo-6)eo)-(6-(eo-6)eo))(r) 

(e 2 (a, 6))(r, «) = c(r, ,) - e (r, ,)(e ■ c) - ei (r, ,)( ei ■ c))(r) (g?) 

V((c- (e • c)e - (ei • c)e x ) • (c- (e • c)e - (ei • c)ei))(r) 

d(r,s) - e (r,s)(e • d) - e 1 (r,s)(e 1 • d) - e 2 (r,s)(e 2 • d))(r) 
(e 3 (a, o))(r, s) = ( 
a/ ((d - (e • d)e - (ei • d)ej - (e 2 • d)e 2 ) • (d - (e • d)e - (ei • d)e x - (e 2 • d)e 2 ))(r) 

We will also define a basis such that /o through / 3 are unit vectors in the directions 
A through D: 



fo{a){p,q) = 7= fi{b){p,q) 



/ 2 (c)(p, 9 ) ^ ^4=, /.(d)(p,ff) - (89) 
V(c-c)(p) y/(d-d){p) 

From this point on we can pursue two different models, which are different from each 
other in terms of the extend to which we apply the context of Ref |T] . In both models we do 
have to use what we refer in the current paper as "part 1" of Ref pQ that deals with literal 
definition of Grassmann numbers since we need that in order to analyze the fluctuations of 
the latter. As was mentioned earlier, even if we choose the integration approach, we would 
still have to analyze the fluctuations of fields in the interior of Alexandroff sets in order 
not to select an Alexandroff set that lies near the lightcone which might result in singular 
behavior due to nonlinearity of fields. However, while we do have to use what we refer in 
current paper as "part 1" of Ref |T] , we don't have to use parts 2 and 3 of Ref [I] . We might 
instead do the same thing as is done in standard quantum field theory: view ifii as degrees 
of freedom independent of vierbines or any other vector fields, and simply couple them to 
vierbines. Thus, we have two options: 

OPTION 1: Use only what current paper calls "part 1" of Ref [1] OPTION 2: Use all 
three parts of Ref [1J 

Both points of view have their own advantages and disadvantages. 

Advantages of options 1: 

a) The measure A involved in option 2 is a lot more complicated than the measure £ 
that will be used in option 1 because of the need to take "derivative" of fermionic degrees 
of freedom with respect to vectors that define them. 
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b) Since in order to talk about fermions we are forced to introduce vierbines in one way 
or the other, it is tempting to use them for a different purpose as well: to adress the question 
as to why causal set manifoldlike. As will be explained in more details in the next section, in 
order for holonomies related to vierbines to be of use in that respect, we would need to use 
random distribution, i.e. measure has to be constant. This is not possible to do in option 
2 since the Grassmann measure A used in the latter takes into account both ip as well as 
holonomies, and A is definitely uneven. In option 1, on the other hand, we can still do that 
since there holonomies and vector fields are treated as separate. 

Advantages of option 2: 

a) If vierbines are treated as separate from the spinor field, they feel like they were put 
in by hand. This make it seem like a certain direction is "better" than others since spin up 
and spin down are taken with respect to that one directoin. The "better" direction is, of 
course, given by vierbines. 

From philosophical point of view, option 2 is clearly more advisable. However, from 
computational point of view option 1 is better since, as you will find, the measure A involved 
in option 2 is a lot more complicated than the measure £ that will be used in option 1 because 
of the need to take "derivative" of fermionic degrees of freedom with respect to vectors that 
define them. Since the ultimate way of testing causal set theories is numerical methods, the 
computational aspect should not be ignored. However, since part of my motivation of doing 
causal sets on the first place is philosophy, the philosophical aspect should not be ignored 
either. For this reason, I resort of doing two different models and separately doing option 1 
and option 2. 



In Option 1, Dirac field ip will be interpretted as 4 independent complex scalar fields. 
Since holonomies a, b, c and d determines vierbines, I will say that these holonomies will be 
interacting with the four components of spinor field I have just introduced. 

The derivative of above defined spinor field is given by 



The gauge field with which it interacts (such as for example electromagnetic field) is 
given by a holonomy v(r, s) which means that in a frame given by our vierbines, 



Option 1 




(90) 



V m = f(e m ,v,a(p,q)) 



(91) 



The interaction of fermionic and gauge field is given by 
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^fy^ = ^^•( 7 V) l i/(e m ,t;,a(p,g)) (92) 
This means that the Lagrangian is given by 



C(a,b,c,d)(p) = ('j°Y)ij^iip)g(em{a,b,c,d),il) j )(p) 

+^i(x P , Xa, a, b, c, d)(p)ip j (p)(>y >y m cr ab ) ij (h(e a (a, b, c, d), e b (a, b, c, d), e b (a, b, c, d))(p) 
-h(e m (a, b, c, d), e a (a, b, c, d),e b (a, b, c, d))(p) + h(e a (a, b, c, d),e b (a, b, c, d),e m (a, b, c, d))(p)) 
+^(p)^h°l%f(e m (a, b, c, d),v)(p) (93) 

Finally, the amplitude is given by 

z = J [da][db][dc}[dd\[dMdMdH[dH {(MfaMMfaMMfoMMfo)) * 

x^(i?e(^))e(/m(^3))e(^e(^ 4 ))e(/m(^))(f 1 A ... A r 8 ) ■ e lS ^ (94) 



Option 2 



In this option, ip will be a function of vector fields, which in causal set case means 
holonomies. We simply rewrite the definition of ip in terms of vector fields given in [1] and 
repeated in equations (14) and (15) and rewrite it in the language of causal sets. This would 
give us the following: 



*Pi(Xp,Xa,a,b,c,d)(p) (95) 
= ( exp { - {(He-'ia, b, c, d)(p)f(a, b, c, d)(p)))^} ).. ( Xp {p)8{ + X a(p)Si) ■ 

Likewise, the measure on ip space is nothing but a copy of equation (16) in causal set 
language: 



Kx P , Xa, a, b, c, d) = £[Re(ip(x P , Xa, a, b, c, d))]£[Im(ip(xp, Xa, a, b, c, d))] x 
xlim^oe^'^x'p, x' a , a', b' , c', d! \i/j(x' p , x' a , a', b', d, d') 
-ip(Xp,Xa,a,b,c,d)\ < e A Va[(e^(a', b', c', d')) - e£(a,b,c,d)) 



(96) 



(e afl (d,b',c',d') 



^afi 



[a,b,c,d))] < e 2 | 
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where \x is a usual Eucledian measure on IR 18 

The Lagrangian density is given by rewriting equation (93) while inserting the definition 

of ij) : 



C(a,b,c,d)(p) = (7 ^) ij ip*(xp,Xa,a,b,c,d)(p)g(e m (a,b,c,d),i) j (xp,Xa,a,b,c,dj)(p) 
, Xa, Q>i b, c, 

p-, Xa-, a , b, c, d)(p)(7 7 m cr a )ij(h(e a (a, b, c, d), e^a, b, c, d), e&(a, 6, c, d))(p) 
-h(e m (a, b, c, d), e a (a, b, c, d),e b (a, b, c, d))(p) + h(e a (a, b, c, d),e b (a, b, c, d),e m (a, b, c, d))(p)) 
+4 ! i(Xp, Xa, a, b, c, d)(p)ip j (xp, Xa, a, c, d)(7°7 M )ii/(e m (a, 6, c, d), v)(p) (97) 

Finally, the amplitude is 

Z = J [da] [db] [dc] [dd] [d X p] [d X a] \{X P , Xa, A, B, C, D)X{ Xp , Xa, A, B, C, D) 

X A(Xp, Xa, A, S, C, ,D)A(X P , Xa, A, 5, C, D)(fi A ... A f 8 ) • e ^(XP,X-^B,C I iJ)) (9g) 

where S(ijj) = f £(ifj)d d x is defined in the usual way, where A product is used for 
multiplication. 



4. Discussion 

This paper completted the main steps of introducing matter fields into causal set theory. 
Scalar fileds were introduced in Ref [3] and Ref [2]. Gauge fields were introduced in Ref [1] 
and Ref [2] and gravitational field was introduced in Ref [3]. The goal of this paper is to 
similarly introduce fermionic fields into causal set theory. 

However, introducing fermionic fields into causal set theory does more than to simply 
complete a list of fields. In particular, there is an important philosophical result of what we 
mean by fermions in general: namely that fermions can be interpreted as local frame defined 
by non-orthonormal vector fields. While similar thing was implied by [lj as well, the lack 
of aforegiven manifold structure in this paper makes this point stronger. In fact, fermions 
can be viewed as agents that cause causal set to be manifoldlike, which is an interesting 
observation since the issue of manifoldlike-ness of causal set is the unresolved one by this 
point. However, one has to be a little more careful. In principle, for any partial order it is 
possible to adjust a set of vierbines in such a way that 6q(p, q) — e\ (p, q) — e^{p, q) — e^(p, q) 
is positive if and only if p and q are causally related. The key thing here is the fact that 
while such choices do exist, if causal set is non-manifoldlike, they are "rare". Thus, since 
integration over holonomies is restricted to the ones that satisfy above condition, the integral 
is small since the range of integration is small. This would imply that manifoldlike causal 
sets are "selected out" as more probable ones. 

The above argument implicitly relies on the fact that measure is uniform in the space 
of holonomies. This would not be true for the "option 2" approach (see above) since in such 
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an approach vierbeins, being "mixed" with spinor fields, are being subject to measure that 
has both positive and negative values. Furthermore, the measure is defined really based on 
the measure of a range of the map from holonomy space to ip space. Thus, if "a lot" of 
holonomies are mapped into the same values of ip, the fact that there are a lot of them won't 
make measure any larger. This means that the fact that in case of non-manifoldlike causal set 
only "few" holonomies would satisfy the restrictions does not imply that the corresponding 
measure is small and thus does not imply the smallness of the integral. This means that if 
we are to use fermions in order to explain the manifoldlikeness of causal set, we are to stick 
to option 1. 

However, option 1 gives us different questions. Since in option 1 fermions and vierbines 
are separate, vierbines are no longer fields. But a vierbine is precisely what tells us what 
direction is the one with respect to which spin might be up or down. So if vierbines are not 
elds, then the question is what does make that direction better? One answer is to go back 
to more standard view in which they are still viewed as fields, but instead of being part of 
fermionic field they are square roots of gravitational field. Incidentally, as already noticed 
by others, this would have an advantage that instead of taking square root of determinant of 
a metric, we would have determinant of vierbine without square root, which might possibly 
adress renormalization problem. 

There are two problems with this approach, however. First of all, if we take this view 
seriously, then we would have to abandon causal structure altogether, as opposed to emposing 
a constraint on consistency of vierbine related holonomies with causal structure. In this case 
we would be forced to adress the question of quantum gravity: how can gravitational field be 
both a stage and an actor? Of course, such question has to be adressed anyway, which is the 
ultimate point of causal set theory. But until it is adressed, we would like to be able to be able 
to describe fermionic quantum field theory on a fixed gravitational background. It might still 
be useful from philosophical point of view to have fixed causal set background as opposed to 
fixed manifold background since I would like to claim that fermions are ultimately responsible 
for causal structure. In this case, the fact that vierbines vary while causal structure is fixed 
might be an argument against viewing vierbines as part of gravitational field, especially in 
light of the fact that option 2 offers an alternative. Although still, in option 1 this can be 
avoided by saying that vierbines do define gravitational field, and they simply vary within 
a certain constraint that el(p, q) — e\{p, q) — el(p, q) — e|(p, q) either stays positive or stays 
negative for any pair of points. 

Another problem with option 1 is that if vierbines are to be viewed as fields, they would 
contradict any other field that we know in a sense that their variation is severely restricted by 
orthogonality and norm 1. This can be adressed through Lagrange multipliers, but viewing 
Lagrange multipliers as fields would only pose more problems. This brings us to the Gramm 
Schmidt process of part 2. This allows us to define vierbines not as fields themselves but 
rather as functions of other fields, and behavior of the latter is not constraint. But in this 
case we would have to worry about the measure due to the fact that we would have to 
take "derivative" of vierbines with respect to original vectors. In this case, one can argue 
that since we have to introduce measure anyway, we might as well introduce the kind of 
measure that is usefull for Grassmann integral, which would bring us right back to option 
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2. Furthermore, since Gramm Schmidt process is not the only way of selecting orthonormal 
basis, one can argue that it is "arbitrary" anyway, and if it is, why can't we "arbitrarily" use 
the remaining degrees of freedom for the definition of fermionic field as is done in option 2? 
Of course, as was said earlier, the price to pay for option 2 is that due to the exotic measure 
we can no longer use fermions to explain the manifoldlike structure of causal set. But one 
can claim that lack of manifoldlike structure is a problem that has always been there, so we 
haven't created any new problems; we simply didn't adress as many problems as we hoped. 
Nevertheless, option 2 DOES define fermions as local frames; it simply doesn't adress the 
question of manifoldlikeness, which means, for example, smooth transition between these 
frame. 
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